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The reflectivity R{ijj) of ab-oriented }Ag\-x^^x{^i-yCy)2 single crystals has been measured by 
means of infrared microspectroscopy for 1300 < ui < 17000 cm~^. An increase with doping of the 
scattering rates in the tt and a bands is observed, being more pronounced in the C doped crystals. 
The CT-band plasma frequency also changes with doping due to the electron doping, while the tt- 
band one is almost unchanged. Moreover, a cr — > ct interband excitation, predicted by theory, is 
observed a.t ujib — 0.47 eV in the undoped sample, and shifts to lower energies with doping. By 
performing theoretical calculation of the doping dependence uib, the experimental observations can 
be explained with the increase with electron doping of the Fermi energy of the holes in the cr-band. 
On the other hand, the a band density of states seems not to change substantially. This points 
towards a Tc reduction driven mainly by disorder, at least for the doping level studied here. The 
superconducting state has been also probed by infrared synchrotron radiation for 30 < < 150 
cm~^ in one pure and one C-doped sample. In the undoped sample (Tc — 38.5 K) a signature of the 
TT-gap only is observed. At y = 0.08 {Tc — 31.9 K), the presence of the contribution of the cr-gap 
indicates dirty-limit superconductivity in both bands. 

PACS numbers: 74.70.Ad, 74.25.Gz, 74.62.Dh 
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I. INTRODUCTION 



Since the discovery of superconductivity in MgB2 
(Tc ~40K)i the electronic properties of this system 
have been intensively studied from a theoretical and 
experimental point of view. The band structure of 
MgB2 is characterized by two distinct electronic bands: 
the quasi two-dimensional a band, formed by the hy- 
bridized spa^Py B orbitals and consisting of two hole- 
like sheets, and the three-dimensional tt band, made 
of Pz orbitals and consisting of two hole-like and one 
electron-like honeycombs iSiSiiiSiSii The disparity between 
TT- and a- band suppresses impurity interband scatter- 
ing giving rise to the most intriguing feature of the su- 
perconductor MgB2, that is multigap superconductiv- 
ity. The (T bands holes are strongly coupled to the in 
plane boron mode E2c,ri2i2ii£ and originate a large su- 
perconducting gap (Act), whereas a small one opens on 
TT-bands (A^),iiJ^J^iil2^J£iiLi8 xhe electronic proper- 
ties of MgB2 can be modified by chemical substitution. 
The two elements which substitute most readily in MgB2 
are Al for Mg and C for B (an exhaustive collection of 
literature data is reported in Ref. [TqI) . Atomic substi- 
tutions change the electronic properties of MgB2 both 
by increasing the impurity scattering (interband and in- 
traband) and by changing the electron density. Indeed, 
C and Al have one more electron than B and Mg re- 
spectively, and therefore it is expected that electrons are 



doped into the system. The most evident effect of Al and 
C doping is the decrease of Tc [see for example Refs. l20l 
andl2ll and references therein]. 

Two changes in the electronic properties with doping 
can be considered. The first change is the shift of the 
Fermi level. Although electron doping is expected from 
the valence consideration, the experimentally observed 
evidence is unclear. The change of the electron den- 
sity was tentatively estimated by means of Hall effect 
measurements but a quantitative and selective analy- 
sis of Hall coefficient is difficult in a multi-band system. 
In MgB2 the band filling due to electron doping should 
reduce the carrier density in the hole like c-band, but 
experimentally it is not yet clear if this effect leads to 
a decrease of the density of states at the Fermi level 
Na{0)- The decrease of N„{0), which enters electron- 
phonon coupling constant A, would lead to a reduction of 
Tc and of the energy gaps A^, A^. The second change is 
the increase of carrier (inter-band and intra-band) scat- 
tering rate. The latter effect was estimated by several 
experiments, as, for example, measurements of resistiv- 
ity and upper critical field iSLSSiSiSiS^ These effects are 
much stronger in C than in A substituted systems. In 
particular, a quantitative estimate of the change of the 
intra-band scattering rate in the tt and a band, selec- 
tively, has been given by far infrared measurement of the 
reflectivity above and below Tc in Al doped and neutron 
irradiated polycrystalline samplesi2& Evidence for the in- 
crease of the inter-band scattering rate was given by mea- 
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suring the merging of the two gap with C substitutioniai 
Indeed, while intra-band impurity scattering is expected 
to reduce Tc by disorder-induced pair breaking, a small 
inter-band impurity scattering may lead to a mix of a and 
TT Cooper pairs, averaging the order parameters (Acr — 
7meV and A^r — 2meV) and reducing Tc down to the 
isotropic value 

Infrared (IR) spectroscopy, when performed on ori- 
ented samples, is a powerful tool to measure both the 
carrier density and the scattering rate selectively in the 
two bands. Furthermore, IR spectroscopy is also a bulk 
probe of the superconducting state, which can provide 
the gap values, as demonstrated recently in polycrys- 
talhne MgB2.2& Although more detailed information can 
be obtained in single crystals than in polycrystalline sam- 
ples, however the small size of available MgB2 single 
crystals (typically 0.3x0.3 mm^) has been up to now a 
strong limit for the use of the IR spectroscopy. Indeed, 
this technique is challenging on such samples, because of 
low signal intensity, further reduced by diffraction effects. 
For this reasons, until now, an extensive IR spectro- 
scopic work at low temperature has been carried out on 
polycrystals^SiSi^Si^i and oriented filmsSiS^ only. How- 
ever, for a deeper insight in the fundamental properties 
of MgB2, oriented single crystal are needed, since in films 
and polycrystals both the intergrain effects and the in- 
crease in the scattering rate of the carriers are present 
and may change the electronic properties qualitatively. 
The infrared response of MgB2 single crystals in the su- 
perconducting state was first investigated by Perucchi 
et alm^ as a function of temperature T and magnetic 
field. To increase the signal in the far-infrared region 
(25 < uj < 50 cm~^), where conventional sources provide 
poor photon flux, Perucchi et al. built a mosaic of several 
crystallites. Although clear effects of the superconduct- 
ing transition could be observed in Ref. [3 the reflectiv- 
ity R{u)) was lower than 70% for 50 < < 6000 cm~^. 
More recently, Guritanu et alm^ investigated the optical 
response of single crystallites of MgB2 in the mid-infrared 
and visible range at room T. They found R{uj) > 90% 
below 6000 cm~^ with field parallel to the a5-plane. This 
value of R{uj) can be reconciliated with band calcula- 
tion predictions^SiS and dc-conductivity {(Jdc = p~^) 
measurementS)2ii2^ indicating that the use of single crys- 
tallites seems to be preferred for an extensive study of 
the electronic properties by infrared spectroscopy. 

In the present work we measure R{oj) in MgB2 single 
crystal at room temperature, and investigate how it is 
modified by C and Al doping. R{uj) in the pure sam- 
ple is in agreement with the results of Ref. |33- Pure 
and doped crystals show a R{uj) metallic, with a pseudo 
plasma-edge at around 2 eV slightly decreasing on dop- 
ing. With increasing C or Al content, the expected in- 
crease of the scattering rates in the tt (F^) and a (F^) 
bands is observed. The plasma frequency of the a band is 
affected by electron doping too, while that of the tt band 
is almost constant. Furthermore, an absorption band at 
~ 0.47 eV is found in the pure sample, which becomes 



less evident and finally disappears as doping proceeds. 
The latter feature, not observed in previous infrared ex- 
periments on pure MgB2 single crystals^ is explained 
here in terms oi a ^ a interband electronic transition 
by direct comparison with band structure calculations. 
The calculated and observed redshift of the cr-interband 
transition with C- and Al- doping allows us to provide 
an estimate of the corresponding Fermi level shift. 

Since the a band in the pure compound fulfills the 
clean- limit condition (Fcr < 2ACT)rithe superconducting 
transition is not expected to strongly modify the infrared 
absorption around huj ~ 2Acr. This is not the case for the 
dirty-limit 7r-band. Indeed previous infrared experiments 
in the superconducting state^^^^i found an increase of 
R{lu) around fku ~ 2A-^ only. However, the increase of 
Fcr is expected to make the effect of the large gap Ao- 
observable in the far infrared spectrum, as demonstrated 
in Ref. ^3 on ^ neutron-irradiated polycrystalline sample. 
In the present work, we show the different far-infrared 
response below Tc of a pure and a C-doped single crystal, 
which can be ascribed to the transition towards dirty 
superconductivity with C-doping in MgB2. 



II. EXPERIMENTAL 

We selected a series of C-doped and a series of Al- 
doped crystallites (see Table IJ). The a6-oriented sur- 
face is not larger than 300x300 fj,m in doped crystals, 
and it often presents cracks and holes when inspected 
at the optical microscope (see Fig^, thus making a 
normal-incidence reflectivity measurement very challeng- 
ing. However, we have selected clean surfaces for reflec- 
tivity measurements in the mid-infrared range by means 
of an infrared microscope, where the beam coming from 
the interferometer is focused in a micrometric spot-size 
by means of a beam condenser. However, infrared mi- 
croscopy is limited on the low-frequency side around 200 
cm~^ (radiation wavelength A of 50 fj,m) by diffraction 
effects. The study of the superconducting state requires 
photon energies of the order of the gap or smaller, i.e. 
uj < 100 cm~^ (< 12 meV). Infrared Synchrotron Ra- 
diation (IRSR) with a conventional normal-incidence re- 
flectivity setup was used in the present work to provide 
a photon flux density high enough to get a signal in the 
far-infrared from a 300x300 /im single crystal inside a 
cryostat. 

For the synthesis of Al and C substituted MgB2 crys- 
tals a high-pressure growth method has been applied^ 
A detailed description of the technique can be found in 
Refs. l23andl2llfor Al and C doped system, respectively. 
The Al and C content of each sample studied is listed in 
Table ^ together with superconducting Tc determined by 
magnetization measurements. 

We used an Hyperion-2000 infrared microscope con- 
nected to a Bruker IF66 interferometer to select clear 
apertures of 40x40 fim on the surface of each sample. 
Gold was evaporated on a portion of the sample surface 
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FIG. 1: Images of three of the samples studied. The large tick intervals correspond to 25 fim. a. Pure MgB2 b. Al-doped 
c. C-doped with evaporated gold (clear region). The white crosses indicate the different points investigated to check possible 
inhomogeneity on the sample surface. The square in the center of the picture is the microscope clear aperture. 



TABLE I: List of sinf 


;le crystals Mgi_ 


xA\x{'Bi^yCy)2 examined in 


the present work, with the 


corresponding code name, C (y) 


and Al (x) content, and superconduct 


ing critical temperature Tc. 






batch 


codc 


X 


y 


Tc (K) 


AN406/1 


P 








38.5 


AN284/9 


C5.3 





0.053 


34.75 


AN286/6 


C8.3 





0.083 


31.9 


AN277/5 


C9.7 





0.097 


29.0 


AN258/2 


C11.4 





0.114 


23.5 


AN262/19 


A7.2 


0.072 





33.05 


AN273/1 


All. 8 


0.118 





30.5 



to get an in-situ reference signal (see Fig. QJ. The fi- 
nal R{uj) is obtained by multiplying for the reflectivity 
of gold. A 40x40 fim clear aperture is large enough that 
there is no need to use the IRSR source, so that a conven- 
tional thermal source and a tungsten filament lamp were 
employed. The reflectivity R{u!) in Fig.|2and|21was mea- 
sured at room temperature for 1300 < ui < 17000 cm^^ 
with a 15x beam condenser. The low-frequency cut-on is 
due to the microscope detector, while the high-frequency 
cut-off, close to the gold screened plasma frequency (~ 
2 eV), is mainly due to the gold fllm deposited on the 
microscope mirrors and to errors in the gold reflectivity 
correction procedure. The response of MgB2 single crys- 
tals in the superconducting state has been investigated 
in the present work by normal-incidence reflectivity ratio 
measurements by FT-IR in the far-infrared range, where 
the photon energy is of the order of the superconducting 
gap {hw < 10 meV). 

The measurements at low-T and low-w were performed 
at the infrared beamline IRIS at the BESSY-II storage 
ring in Berlin, Germany. The beam collected at the end 
of a bending magnet from vertical angle of 40 mrad is sent 
to a Bruker IF66 interferometer equipped with a liquid- 
He cryostat and a normal-incidence reflectivity setup. A 
cut-on optical fllter at 100 cm~^ was put in front of a 4.2 
K bolometer. The intensity of the IRSR source is about 



one order of magnitude larger than that of a Hg-arc lamp, 
depending on the current inside the storage ring (typical 
values are 150 250 mA). However, the signal-to-noise 
ratio of the far-infrared data from a 300x300 /im sam- 
ple is strongly limited by several experimental problems, 
which we now discuss. Information on the gap value 
of MgB2 has been obtained from FT-IR measurements 
in far-infrared on pelletS)S& filmsj^^ and single crystal 
mosaics^ However, we decided to measure the intensity 
reflected by one single crystal, not a mosaic or a film, in 
order to be able to compare the data with quantitative 
models of the electrodynamic in-plane response, based 
on the Bardeen-Cooper-Schrieffer (BCS) theory, which 
includes the plasma frequencies, the scattering rates and 
the superconducting gaps of both the tt and the a band 
(see section IV). The disadvantage is that the available 
surface area for the rcfiection (300x300 ^m^) is 50 to 1000 
times smaller than that, for example, of Refs. 12^13313^ 
The loss of signal intensity could in principle be recov- 
ered by the use of infrared synchrotron radiation (IRSR) , 
which provides a much higher brilliance if compared with 
blackbody sources. However, two main problems arise for 
u! 0: diffraction effects and increase in the IRSR beam 
size. In the far-infrared range (20 < oj < 100 cm~^, 
2.5 < E = fku < 12.5 meV) the radiation wavelength 
A is 0.1 < A < 0.5 mm, of the order of the sample size 
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d ~ 0.3 mm. Therefore, we are deeply into the diffraction 
regime, where a large portion of the incoming radiation 
is not reflected at specular angle, but rather diffused in a 
much larger solid angle. The effect of diffraction on the 
measured signal is difficult to be exactly determined, as 
it depends on the details of the experimental setup geom- 
etry. We have however estimated, through a scries of test 
measurements in our setup, that diffraction effects should 
decrease the signal intensity by one order of magnitude 
in the entire far-infrared range4£ The second problem is 
that the IRSR beam size approximately increases as X^^^ 
(Schwinger law). Since the spot size at the sample posi- 
tion in the visible range is about 500 /xm in diameter, the 
spot at 30 cm^^ should be around 5 mm in diameter4i 
As a consequence of this, only a minor part of the IRSR 
photons hit the sample. 

To eliminate the diffraction effects, we eliminated any 
uncertainty on the optical alignment due to cryostat 
strain or to IRSR variations by directly measuring the 
reflected intensity ratio (which in principle equals the re- 
flectivity ratio) Rs/Rn = R{uJ,T = 4:.2K)/ R{lu,T = 
45K) during several experimental runs without heating 
above T = 50 K. In this way, the diffraction pattern is 
washed out by dividing the spectra taken with the same 
setup geometry. To avoid the contribution from the ra- 
diation not hitting the sample, we mounted the crystals 
on the top of a cone-shaped sample holder. To check 
the contribution from diffuse radiation, we performed a 
background intensity test by measuring the intensity re- 
flected by the sample holder without the sample. The 
background signal is smaller by a factor 20 and docs not 
change for cryostat temperatures T < 50 K, and there- 
fore could be easily subtracted. In conclusion, we could 
obtain reliable data for 30 < a; < 100 cm~^ by averag- 
ing a total of 20000 interferometer scans, with a spectral 
resolution of 2 cm~^ and an uncertainty of ±0.5% on 
the reflectivity ratio. We repeated the procedure several 
times for one pure and one C-doped sample (y = 0.083). 

III. NORMAL STATE 

The reflectivity i?(a;) at room temperature of the two 
single crystal series is shown in Fig. [5| and 13 in the range 
of the present measurements (0.2 < w < 2 cV). For all 
crystals, R{w) is a typical free carrier response in a metal, 
with a pseudo plasma-edge (i.e. the frequency where 
R{uj) almost reaches zero value) around 2 eV and almost 
independent on doping. The pseudo plasma-edge roughly 
corresponds to the screened plasma- frequency u)p, which 
is determined by the free carrier density in both the a 
and TT bands by means of the relation: 

= (1) 

where is the high frequency dielectric constant deter- 
mined by the electronic interband transitions and uj^ , lUj^ 



E(eV) 

0.0 0.5 1.0 1.5 2.0 



^1 1 1 
- ^ p 


1 






~ ■■ — - — ^..coX. \, 


- 




w 






- Mg(Bi_yCp2 




single crystals 




~ T=300 K 




1 1 



4 8 12 16 

(0(10^ cm"b 

FIG. 2: Reflectivity at 300K of C-doped crystals with field 
parallel to the ab plane. The curves are shifted for clarity. 



are the plasma frequencies respectively for the a and tt 
bands. We note that the level of the reflectivity in the 
mid-IR (w < 1 eV) decreases with increasing both Al and 
C doping, indicating the expected increase in the scatter- 
ing rates T^, F^ (please note that the curves in Fig.|21and 
Oare shifted for clarity). In Fig. [3 and El a kink at ~ 0.4 
eV in R{uj) is evident in the pure sample. This feature, 
not observed in previous infrared experiments on pure 
MgB2 single crystals, becomes less prominent and finally 
disappears as doping proceeds. 

Since there are several open questions about the qual- 
ity of optical data on MgB2 crystals, it is worth to 
compare in Fig. our reflectivit y sp ectrum on pure 
MgB2 with those reported in Refs. l34l and [3^ The lu- 
dependence of the free carrier contribution is remarkably 
similar in all three cases. A small discrepancy is seen 
below ~ 0.5 eV, where we observe a change of slope, 
not so clearly seen elsewhere. One may ask whether the 
present observation can be ascribed to surface contami- 
nation elfects. Although the reason of this discrepancy is 
not clear at the moment, however, the reflectivity level in 
the present work is very close to that of Ref.|3^ where the 
surface contamination is minimal. Another discrepancy 
between our data and those of Ref. [s^ is seen above 2 eV 
and probably due to an artifact of the gold reflectivity 
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FIG. 3: a. Reflectivity at 300 K of Al-doped crystals with fleld 
paraUel to the ab plane. The curves are shifted for clarity, b. 
Comparison of Riuo) of pure MgB2 single crystals from present 
work and from two previous works i—*— 
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FIG. 4: Reflectivity of pure, C-, and Al-doped MgB2 crystals 
with field parallel to the ab plane in the low frequency range. 
The thin lines are the result of the fitting procedure described 
in the text. 
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correction described in section I. 



The analysis of the optical spectra allows for a reliable 
determination of the plasma frequencies and the scatter- 
ing rates. Considering that it seems very hard to ex- 
tract quantitative values for and uj„ from the mea- 
sured screened plasma frequency, since the exact varia- 
tion of with doping is unknown, we employed a Drudc- 
Lorcntz fitting procedure, which was shown to work well 
in the case of pure MgB2i^ The Drudc model can be 
adapted to the case of MgB2 by introducing two sepa- 
rate Drude terms for the two bands, plus a suitable £00 
which we take as doping independent and set to the value 
of 11.9 found by ellipsometric measurements in Ref. Issl 
A Lorentzian oscillator centered at lomir = 0.47 eV was 
necessary to reproduce the change of slope seen in the 
reflectivity. The final fitting formulas arc therefore 



1 - V~e{uj) 



1 + 



(3) 



where Smir,^ mir. a-re the strength and the width of 
the Lorentzian oscillator respectively. In Fig. ^ we show 
the result of the fitting procedure together with the re- 
flectivity data, which are well reproduced by the above 
formula. The discussion about the evolution with doping 
of the resulting fitting parameters is detailed in the next 
section. 

The optical conductivity a{ijj) was obtained from R{lij) 
through standard Kramers-Kronig transformations. The 
reflectivity was extrapolated to zero frequency by using 
the fitting curves in Fig. 0] below 1300 cm~^. To extrap- 
olate at high frequencies wc used the data from Ref. ISSi 
showed in Fig. ^p. In Fig. |S1 where <j{lo) is plotted for 
representative doping values of C (panel a) and Al (panel 
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FIG. 5: Optical conductivity of selected single crystals. The 
arrows indicate the value of lomir as obtained from fitting 
the reflectivity data (see text). 



b), a typical metallic Drudc-like behavior is clearly evi- 
dent at low frequencies for all samples. The extrapolation 
of a(ui) to zero frequency strongly decreases with doping 
in agreement with dc transport measurements on similar 
samplesiSi At mid-lR frequencies, between 2000 and 5000 
cm~^, the deviation from a monotonous behavior, more 
pronounced in the pure sample, indicates the presence 
of a broad absorption band superimposed to the Drude 
peak. The black arrows indicate the frequency ujmie. of 
the Lorentzian oscillator at 0.47 eV introduced in the fit 
to the reflectivity data. The overall agreement between 
^MiR obtained from the fit to R{u)) and the bump in the 
optical conductivity indicates that the change of slope in 
R{ijj) corresponds to the mid-IR absorption band. The 
latter feature clearly shifts and broadens as the doping 
(Al and C) increases. In the discussion section of this 
paper we will argue that this peak corresponds to an in- 
terband electronic transition between two a bands. 



A. Discussion 

The main effects of Al and C substitution are expected 
to be an increase of the scattering rate and a decrease of 
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FIG. 6: Scattering rates Fo-jF^ as obtained from the fit of 
Eq. 13 to the refiectivity data, for C (panel a) and Al (panel 
b) doped crystals. The dotted lines are guides to the eyes. 



the plasma frequencies in the a and tt bands. Looking 
at the FiglHI we notice a clear evolution with doping of 
the scattering rate in the two bands, while the plasma 
frequency seems to vary only in the a-band. Let's first 
discuss the evolution of the scattering rates with doping. 

In Fig. Elwe show F^ and F^r in the two series of Al and 
C doped samples as a function of dopant concentration, 
that is the amount of substituted Mg or B atoms. 

In the pure sample Fo- = 37 meV and F^r = 87 meV 
in good agreement with the values found in Ref. IbM 
This finding (Fq- <C F^) is expected, since the cr-band 
is considered in a clean conduction regime, whereas the 
TT-band in a dirty oneiiiS^ With increasing Al and C dop- 
ing the scattering rate increases in all samples of the two 
series, but with a different trend. In the C doped sam- 
ples Fo- and F^ both sharply increase towards a common 
value at high doping level. This means that C substitu- 
tion is much effective in disordering both the bands, but 
more strongly the a-band driving it towards a less clean 
conduction regime. This is consistent with the strong 
increase of the upper critical field measured in C-doped 
single crystals iSiiS^ However, a different scenario could 
also be considered. In the present fitting procedure, we 
can distinguish between the two bands only if they have 
a very different scattering rate. Therefore, we cannot ex- 
clude that for high C- doping levels, the conduction takes 
place in one band only. As for the Al doped samples, the 
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FIG. 7: a. Plasma frequencies of tlie tt (triangles) and a (cir- 
cles) bands and of the total plasma frequencies (squares) as 
a function of electron doping, as obtained from the fit of the 
Eq. I^to the reflectivity data. b. Electron doping dependence 
of the absorption band given by the resonance frequency of 
the Lorentzian oscillator present in Eq. |H] The grey thick 
line is the calculated doping dependence of the a —> cr in- 
terband transition (see text). Dashed line is a guide to the 
eyes. Empty symbols in both panels refer to Al-doped sam- 
ples, filled symbols to the pure and C-doped samples. 



effect of doping on the intraband scattering rates is simi- 
lar in the two bands: they increase substantially only for 
X > 0.07, keeping < ^tt- This result suggests that C 
substitution is more efficient than Al one in disordering 
both the bands. 

In addition to tunc the magnitude of the disorder scat- 
tering processes, Al and C substitutions arc expected 
to introduce also additional charges in MgB2 affecting 
thus the plasma frequency of the system. The measured 
plasma frequencies iOp^a, ^p,Tr, obtained from the fitting 
procedure of the reflectivity data, as well as the "total" 

plasma frequency ujp = 



tot, 



■ iOp .^, are shown in Fig. 

|3i as a function of electron doping for all the samples 
studied (full symbols refer to C doping and empty sym- 
bols to Al doped crystals). 

For the undoped compound we find ujp^^ = 2.6 eV, 
^p.TT = 4.2 eV, yielding ujp = 5.0 eV, much larger than 
previously reported in optical experiments, apart from 
Ref. [S^ The discrepancy between the present and the 
previous measured values of the plasma frequency can 



be ascribed to different sample purity and air contam- 
ination, as also proposed in Ref. l35i In the present 
study, the samples are single crystals with minimal ex- 
position to air, which explains the much larger value. 
However, the present values, in particular iOp^^ = 2.6 eV 
and Up = 5.0 eV, are still smaller than those reported 
in Ref. |3^ {uip ~ 6.3 eV) and those predicted by the 
theory {ujp ~ 7 eVSiLiS). The discrepancy between the 
theoretical and optically measured values of the plasma 
frequencies has been addressed in Ref. [s^ We can note 
however that the discrepancies between the present val- 
ues and the ones reported in Ref. can be explained by 
the difference in the fitting procedure, since we added an 
absorption at 0.47 eV whose spectral weight is compara- 
ble with that of the Drudc terms. Different hypotheses 
can be invoked about the origin of this optical structure. 
We cannot exclude for instance that a minimal exposi- 
tion to air is responsible for shifting part of the spectral 
weight from the Drudc term to finite frequency absorp- 
tion corresponding to charges localized by defects in the 
crystal structure produced by air exposition. In this case, 
part of the spectral weight of the band at 0.47 eV should 
be added to the computation of ujp and a value around 6 
eV could be easily attained. However, as we will see later, 
we propose for the band at ~ 0.47 eV a different physi- 
cal meaning related to a predicted interband transition, 
as suggested by the doping dependence of its frequency 

As for the doping dependence of ujp^i, (v = cr, tt), we 
note the tt band plasma frequency does not change sub- 
stantially within this electron doping range, whereas the 
a band presents a slight decrease. The reduction of 
LOp^rr can be easily explained in terms of reduction of 
the hole charge density in the a bands which, in this 
doping regime, can be safely modelled as almost two- 
dimensional parabolic bands. A simple calculation gives 
thus Wp,CT oc y/Ncr{0)Ep/m*, where Ep is the hole Fermi 
energy of the a bands which is roughly given by the dis- 
tance between the chemical potential and the top edge 
of the (J bands. Theoretical calculations show that the 
a band density of states Na-{0) is essentially doping in- 
dependent in this range of doping^ii so that the de- 
crease of ujp,a- is mainly ruled by Ep. According to 
this analysis, we can thus estimate a reduction of Ep 
of a factor ^ 0.36 from the pure to the highest doped 
compound, which, using Ep ~ 0.6 eV for the undoped 
system^ gives E§,{C11A) ~ 0.22 eV for the C11.4 com- 
pounds, corresponding to s shift of the chemical potential 
= AE"^ ~ 0.38 eV. Note that, since the Fermi en- 
ergy of the TT bands is much larger than the a band one, 
Ep ^ 5 eV, the relative change in the tt band charge den- 
sity is much smaller and practically negligible, in agree- 
ment with Fig. 12^. We have to warn however that some 
caution should be employed when extracting quantitative 
information about ujp^^ and cOp^^ from the Drude-Lorentz 
fit in the high doping region. In this regime indeed the a 
and TT band scattering rates Fg., F^ become very similar, 
making hard to separate in the fit procedure the two dis- 
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tinct a and tt band contributions, whereas of course the 
total value of the plasma frequency ujp is not affected by 
this kind of problem. 

Although a quantitative analysis of the ujp^„, ujp.-K 
would requires some care in the high doping range, a 
clearer evidence of the filling of the a bands as a func- 
tion of the electron doping may be given by the doping 
dependence of the resonance frequency of the Lorcntzian 
oscillator present in Eq. O We have already mentioned 
above that a possible explanation for this structure could 
be related to the interband a ^ a electronic transitions. 
An absorption band around 0.4-0.7 eV has been indeed 
theoretically predicted in the pure sample, although not 
observed so iaxm^i^ However, there could be other dif- 
ferent hypothesis to consider. One already mentioned 
possibility is the presence of an absorption correspond- 
ing to charges localized by defects in the crystal structure 
produced by air exposition. But the high values of both 
reflectivity data and plasma frequency seem to contra- 
dict this hypothesis. Another possible interpretation of 
the MIR band could be given in terms of the frequency- 
dependent electron-phonon scattering rate t~^{uj). In- 
deed, if the conductivity is calculated from the Eliash- 
berg electron-phonon spectral function^^k then a devia- 
tion from the simple Drude model in Eq. |31 is found in 
the MIR range, which can be interpreted as an increase of 
t~^{ll!). However, due to the addition of dopant atoms, 
the scattering should be dominated by impurity scatter- 
ing at all Lu (apart from that in the cr-band of the pure 
sample) and hence the conductivity could be as well de- 
scribed by a frequency-independent t^^ = F in Eq.|21plus 
a MIR band, as for example in the case of cupratcsi^ 
The presence of the MIR band with a similar intensity 
in all samples suggests therefore that it might not be 
fully explained in terms of electron-phonon scattering. A 
last hypothesis could ascribe this MIR absorption band 
to a possible presence of MgO impurity on the sample 
surface. The clear doping dependence of the absorption 
frequency lomir shown in Fig. makes us confident 
in excluding also this last possibility. Indeed, the soft- 
ening of the ujmib. with doping, as shown in Fig. [7|3, 
is expected since the distance between the two hole-like 
quasi-parabolic a bands decreases as electron charges are 
added to the compound. 

To have a more quantitative insight, we compare these 
data with the simple model for the a bands introduced 
in Ref. 0, namely 

ek.i = Co - [kl + kl)/m, - 2t± cos{k^c), (4) 

where i = 1, 2 is the index of the two a bands. Following 
Ref. Qwe take eq 0.6 eV, t±_ ^ 0.092 eV, mi = 0.20TOe, 
TO2 = 0.53TOe- Within this model the a band IB absorp- 
tion results to have a energy window 0.26 eV < w < 1.29 
eV, and a Lorentzian fitting gives u)ib — 0.51 eV, in very 
good agreement with our data. We can also investigate 
the doping dependence assuming a rigid band filling. It 
is easy to show that w/b oc Ep . The results are shown in 
Fig. \7}p, also in an overall agreement with our data. The 



slight bend of our data with respect to the almost lin- 
ear behavior of the theoretical analysis can be probably 
ascribed to non-rigid band effects, which are also known 
to be important in the doping range>i^ Note that the re- 
duction of the measured values of ujib , from the pure to 
the highest doped compound, is ^ 55%, also in a qual- 
itative agreement with the reduction of Ep obtained by 
'^p,a-- We conclude therefore that the assignment of the 
~ 0.45 eV absorption with the IB a transition is in good 
agreement with theoretical predictions, thus giving an ex- 
perimental indication of an upwards shift of the chemical 
potential with the substitution of Al and C and of the cor- 
responding electron doping of the a bands. As discussed 
before, since the a bands are almost two-dimensional and 
the corresponding DOS almost flat in this doping regime, 
the band filling is however unlike to be the only source of 
the Tc suppression, and we think that an important role 
is played by the disorder as well. 



IV. SUPERCONDUCTING STATE 

The reflectivity ratio {Rs/ Rn) is reported for the pure 
and the C-doped sample in Fig. [Hj-a. The pure sample 
shows a decrease below 1 for 40 < w < 80 cm~^ and a 
sharper increase above 1 for a; < 40 cm^^. The C-doped 
sample only shows an increase above 1 in Rs/Rn for 
(jj < 60 cm^^. Above 80 cm~^, there seems to be no effect 
of the superconducting transition on the Rs/Rn within 
uncertainties. Although the data quality does not allow 
an uncostraincd best-fltting procedure, the fact that both 
the present datasets are each taken on one oriented sin- 
gle crystal poses strong constraints. We are therefore in 
the condition to compare the present data from MgB2 
samples with the prediction of different BCS-based mod- 
els for the infrared response of a two-gap superconduc- 
tor. The BCS calculation employed here for the complex 
conductivity a of one band in the superconducting state 
is that proposed by Zimmermann et al.m^ Therein, the 
normal state conductivity is a Drude term defined by the 
plasma frequency LUp and the scattering rate F, while the 
superconducting state conductivity is computed trough 
the input of T/Tc and the gap value A. 

However, in the present experiment the quantity which 
is probed is the reflectivity rather than (7{lu). In 

a s-wave BCS superconductor, for Hui < 2A, the R{uj) 
in the superconducting state is close to 100%, because of 
total radiation screening from the supercurrent flowing 
in a surface sheet, whose thickness is determined by the 
field penetration depth A. The normal state reflectivity 
is then recovered even at T <C Tc for Hlu > 2A, since the 
radiation with energy huj larger than that of the Cooper 
pairs 2A cannot be screened. If one now turns to the 
two-gap case, which is relevant for our work, one may 
ask whether the radiation would be screened or not for 
2A^ < hu) < 2 Act, since the photon energy is now larger 
than the 7r-band Cooper pair energy, but smaller than the 
(T-band Cooper pair energy. Normal electrons, which are 
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FIG. 8: (a). Reflectivity change across the superconducting 
transition in one pure (circles) and one C-doped {y — 0.083, 
triangles). The symbol size corresponds to the estimated stan- 
dard deviation of the spectrum average. (b,c). The data from 
pure (b) and C-doped (c) crystals in a. are shown separately 
and compared to two BCS model curves (see text) . Thick line: 
gap parameters in case 1 (parallel conduction), dashed line: 
gap parameters in case 2 (single effective gap) . (d) The model 
reflectivity curves which generate the Rs/Rn curves in (b,c). 
Thin line: normal state; dashed and thick lines represent the 
superconducting state in the case 1 and 2 respectively. 



unpaired due to radiation-induced Cooper-pair breaking 
in the 7r-band, will coexist with the surface supercur- 
rent generated by the pairs in the cr-band. The question 
whether the radiation with 2At^ < Huj < 2Aa- would be 
screened or not has not yet been addressed theoretically 
to our knowledge. We considered these two possibilities 
by using two different sets of parameters to calculate the 
conductivity in the two bands (case 1 and case 2 de- 
scribed below). We then summed up the contribution 
from the two bands and we determined the reflectivity 



and the Rs / Rn from the complex conductivity with the 
standard Fresnel formulae. 

In case 1 we assumed that the conductivity in the su- 
perconducting state is the parallel of two superconduc- 
tors, each one with a different Wp.^ , and Ai,(p = cr, tt). 
This assumption, which holds for the normal state, was 
already used to analyze the far-infrared data in the su- 
perconducting state of MgB2 pellets and films p^S*^ The 
superconducting state conductivity is then given by sum- 
ming up the contributions from the two bands, each 
weighted by tOp^^. The fact that <C 3.52kBTc must 
be explicitly taken into account in the finite temperature 
calculation. The use of the parallel sum in the supercon- 
ducting state corresponds to the assumption that there 
is no interaction between the two superconducting flu- 
ids when they are exposed to an incoming electromag- 
netic wave. In case 2 we assume that the unpaired tt- 
electrons inhibit the radiation screening for frequencies 
2A^ < hu! < 2 Act, thus allowing for the propagation of 
the radiation in the sample in a surface layer larger than 
the theoretical value of A^. In this case, the reflectivity 
measurements are only sensitive to the smaller gap A^. 
We have therefore input in our calculations one single 
effective optical gap A = Ajr for both bands in case 2. 
We note that the above discussion does not imply that 
there is a single value for the gap or for the penetration 
depth in the two-band superconductor, but rather one 
single critical energy of the order of 2A^ for the ac-field 
penetration depth. 

We have generated two curves for each of the two sam- 
ples, reported in Fig. (Sj-b and -c, by using different gap 
parameters for case 1 and case 2. The input parame- 
ters, reported in Table 2, are entirely determined a pri- 
ori, since the data quality does not allow a reliable fitting 
procedure, c^p ^r and Wp^o- were taken from the room tem- 
perature optical data reported in the present work on 
the same samples, and assumed to be temperature in- 
dependent. The scattering rates are mainly determined 
by impurity scattering at low temperature. Therefore, 
we used the room temperature values for all the cases 
where impurity scattering dominates the electron scat- 
tering also at room temperature, i.e. F^ for the pure 
sample and F^,Fa- for the C-doped sample. According 
to several authorsjSiiSSiS^ in the pure sample the elec- 
tron scattering in the cr-band is dominated by phonons 
at room temperature, so we cannot use the value for Fo- 
from the present work. However, in Ref. |3^ an impurity 
scattering rate value of 12.4 meV was derived and found 
to be in agreement with dc transport measurements; we 
have therefore used this value for the low-temperature 
of the pure sample. Concerning the gap values A^ and 
Act, we have used the values from photoemission mea- 
surements of Ref. for samples with similar doping 
level. T/Tc was 0.1 for the pure sample and 0.13 for the 
C-doped sample. 

The result of the analysis of Fig. |S1 is the following: 
case 1, i.e. the parallel conductivity (dashed line), can 
account for the data from the C-doped sample but not for 
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TABLE II: Input parameters used to generate the model 
curves in Fig. El a nd El All data are from the present work, 
apart " from Ref. l43l and from Ref. Issl All values are in 
meV, where not stated otherwise 
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those on the pure sample. On the other hand, case 2, i.e. 
the single effective gap case (thick line) may explain the 
data from the pure sample. This is not surprising, since 
all the infrared reports of the superconducting gap value 
in undoped MgB2 are in agreement with the value for 
from other techniques i^SiS*^ This low value of the 
optical gap was often explained with the fact that the 
CT-band is a superconductor in the clean limit (F < 2A) 
while the 7r-band is in the dirty-limit (F > 2A). In a 
clean-limit superconductor, the superconducting transi- 
tion would not leave imprints in the electrodynamic re- 
sponse at the gap energy^ and therefore infrared radi- 
ation would only probe the smaller 7r-gap. On the other 
hand, the data on the C-doped sample, with both bands 
in the dirty limit, can be successfully reproduced by a 
parallel sum of two dirty superconductors. As a conse- 
quence, we can confirm that the gap in the clean-limit 
CT-band is not observed in the infrared experiments. On 
the other hand, even a fit curve which includes the clean- 
limit condition can not reproduce the data for the pure 
sample. In the absence of more refined theoretical models 
for the electrodynamic response of a two-gap supercon- 
ductor, we speculate that the lack of CT-gap signature in 
the reflectivity in the pure sample could be due to un- 
paired TT-electrons inhibiting the radiation screening for 
frequencies 2A7r < huj < 2Acr. This is not the case when 
the CT-band is also in the dirty limit, like in the C-doped 
sample, as shown by the higher cutoff frequency in the 
Rs/Rn- We believe that the larger A, and hence the 
larger thickness of the screening supercurrent sheet, may 
play a role in determining the higher cutoff frequency in 
the C-doped sample. Note that, in spite of the higher 
cutoff frequency in the Rs/Rn, the absorptivity 1 — i? of 
the dirty C-doped sample is larger than that of the pure 
sample, as expected (see Fig. (Hj-d). 

The clean-limit nature of the CT-band in the pure sam- 
ple can be made clear by looking, in Fig. IHl at the 
model optical conductivity which was used to generate 
the Rs I Rn curves in Fig. |S| The normal-state conduc- 
tivity in Fig.Elis the sum of two purely Drude terms, and 
therefore cannot be taken as a description of the real fre- 
quency dependence of the optical conductivity of MgB2 
single cystals. However, since the model curves in Fig. 
1^1 are obtained from a-priori values of the input parame- 
ters and generate Rs/ Rn curves in good agreement with 
the experimental data, they can be interestingly used to 
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FIG. 9: (a). Model optical conductivity of the pure sample in 
the normal (thin line) and superconducting (thick line) state 
as determined from case 2 gap values (see text), (b). The 
same as in (a) but for the C-doped sample, thick line obtained 
from case 1 gap values, (c). The spectral weight difference 
between the normal and the superconducting state S (right 
scale) for both the pure sample in (a) and the C-doped sample 
in (b) . The thin dashed lines indicate the asyntotic value and 
90% of the asyntotic value for both samples. On the left 
scale the corresponding value of the squared superconducting 
plasma frequency is reported. S of the pure sample shows 
a clean-limit behavior, since it reaches the 90% line for ~ 
2 Act, while the C-doped sample reaches the 90% line at ~ 
6 Act, hence showing a dirty-limit behavior. 



estimate for both samples the spectral weight loss below 
Tc- This quantity is given by 



(5) 



where cti^„, cti^^ are the real part of the normal and super- 
conducting state conductivity respectively. The quantity 
5(oo), i.e. the value of the integral in Eq.|Slfor lj — > oo, is 
proportional to the superfluid density, as the total spec- 
tral weight lost through the transition is transferred to 
the zero-frequency delta-function representing the con- 
tribution of the superconducting carriers to the optical 
conductivity (Ferrel-Glover-Thinkham sumrule). 

However, in our BCS approach, where the energy scale 



11 



of the superconducting phenomenon is set by the larger 
energy gap Ag., 90% of S{oo) can be recovered by ex- 
tending the integral only up to few units of A^- . as shown 
by the dashed lines in Fig.|51-c. In the model for the pure 
sample (Fig.|51-a) we find that 90% of S{oo) is obtained by 
integrating up to 12 meV ~ 2Aa only. This corrcpsonds 
to the clean-limit case, since most of the normal-state 
spectral weight is located in the frequency region w < 2A 
because of the small scattering rate < 2A, which de- 
termines the width of the Drude peak. On the other 
hand, the model for the C-doped sample shows a broad 
Drude term and therefore a recovery of 90% of S{oo) by 
integration up to 25 meV ~ BA^ (being Acr = 4.2 meV 
in the C-doped sample) which is usually considered the 
upper limit of the integration in Eq.Elfor dirty-limit BCS 
superconductors^^ 

S{oo) can be more readily expressed in terms of a 
superconducting plasma frequency ujps = \J (2/7r)5(oo). 
LOps can be compared to the total plasma frequency tjp 
derived from the fit in Section 2 of the present work to 
give an estimate of the proximity of MgB2 to the ideal 
clean-limit condition, defined by the London model pre- 
diction LOys = at T ~ 0. We obtained lo^s = 1.8 and 
0.9 cV for the pure and C-doped sample respectively (see 
Fig. El-c), to be compared with ld^ = 5.0 and 4.7 eV re- 
spectively. We can conclude trough the analysis of the 
model conductivity curves that the Rs/Rn far- infrared 
data in Fig.|Hlconfirm the clean- limit scenario for the pure 
sample and the the dirty-limit scenario for the C-doped 
sample with y = 0.083. 

V. CONCLUSIONS 

A systematic investigation of the in-plane optical prop- 
erties of Mgi_^Ala;(Bi_.yCy)2 sluglc crystals has been 
curried out as a function of Carbon (x) and Aluminum 
(y) concentration by means of infrared microspectroscopy 
for 1300 < oj < 17000 cm~^. For pure and doped crystals 
the rcfiectivity R{uj) is metallic, with a pseudo plasma- 
edge at around 2 cV slightly decreasing on doping. Us- 
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ing a Drude-Lorentz fit (see Eq. to the experimental 
R{lo), an increase of the scattering rates in the tt (F^r) 
and (T (Fcr) bands is found with increasing Al and C con- 
tent. Whereas the 7r-band plasma frequency results not 
to be affected by electron doping, the a-band one shows 
a substantial decrease. An absorption band at ^ 0.47 eV 
was found in the pure sample, which becomes less evident 
and finally disappears as doping increases. We performed 
band structure calculations which ascribe it to a cr ^ cr 
interband electronic transition. The calculated and ob- 
served redshift of this transition with C- and Al- doping 
allows us to provide an estimate of the corresponding 
Fermi level shift, which is in qualitative agreement with 
the decrease of the cr-band plasma frequency. Since the a 
bands are almost two-dimensional and the corresponding 
DOS almost fiat in the doping region studied here^^ii the 
band filling is unlike to be the only source of the Tc sup- 
pression, and we think that an important role is played 
by the disorder as well. 

The effect of doping in MgB2 has been also probed in 
the superconducting state by using infrared synchrotron 
radiation for 30 < a; < 150 cm~^ in one pure and one 
C-doped single crystal. The far-infrared response below 
Tc in the two samples is substantially different, due to 
the different cr-band conduction regime. Indeed, in the 
undoped sample, with clean cr-band, a signature of the 
TT-gap only is observed. In the C-dopcd one, a contri- 
bution from the cr-gap to the refiectivity ratio appears, 
indicating a transition towards dirty superconductivity 
with C-doping in MgB2. 
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